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1. Introduction 

In this work we present the first lattice QCD evaluation of the N to A axial form factors. A 
determination of these form factors provides an important input for the GO experiment, which is 
under way to measure these form factors at Jefferson Lab The interest in the axial N to A 
transition arises from its purely isovector nature, which probes different physics from what can be 
extracted from the study of strange isoscalar quark currents. Combining results from the electro- 
magnetic N to A transition we evaluate the dominant contribution to the parity violating asymmetry 
as determined by the ratio Cj/C^. This is the analog of the gA/gv rat io extracted from neutron 
/3 -decay. Furthermore we investigate low-energy consequences of chiral symmetry, such as the 
non-diagonal Goldberger-Treiman relation. 

2. Computational aspects 

Given that this is a first lattice computation of the axial transition form factors we test our 
techniques in the quenched theory where we can use a large volume and have small statistical er- 
rors. The large spatial size of the lattice allows to both reach small q 2 values as well as extract 
the q 2 -dependence more accurately having access to more lattice momentum vectors over a given 
range of q 2 . Pion cloud contributions are expected to provide an important ingredient in the de- 
scription of the properties of the nucleon system. In this work the light quark regime is studied 
with pion masses in the range of about (690 — 360) MeV using two degenerate flavors of dynami- 
cal Wilson configurations ^ ||] and in a hybrid scheme which uses MILC configurations generated 
with staggered sea quarks |Q] and domain wall valence quarks that preserve chiral symmetry on 
the lattice. An agreement between the results from these two different lattice fermion formulations 
provides a non-trivial check of lattice artifacts. In particular finite lattice spacing, a, effects are 
different: both the quenched and unquenched Wilson fermions have discretization errors of &(a), 
while both Asqtad and domain wall actions have discretization errors of &{a 2 ). Furthermore do- 
main wall fermions preserve chirality, in contrast to Wilson fermions. The hybrid calculation is 
computationally the most demanding. The light quark domain wall masses are tuned to reproduce 
the mass of the Goldstone pion of the staggered sea. Throughout this work the bare quark masses 
for the domain wall fermions, the size of the fifth dimension and the renormalization factors Za for 
the four-dimensional axial vector current are taken from Ref. [Qj. In all cases we use Wuppertal 
smearing [|]] for the interpolating fields at the source and sink. In the unquenched Wilson case to 
minimize fluctuations [Q] we use HYP smearing [^] on the spatial links entering in the Wuppertal 
smearing of the source and the sink whereas for the hybrid case all gauge links in the fermion 
action are HYP smeared. In Table 1 we give the parameters used in our calculation [^J. The value 
of the lattice spacing is determined from the nucleon mass at the chiral limit for the case of Wil- 
son fermions whereas for the hybrid calculation we take the value determined from heavy quark 
spectroscopy flT0|]. 

3. Methodology 



The calculation of the axial form factors makes use of the same methodology as the one used 
in our lattice study of the electromagnetic N to A transition [fiTl [12|]. The invariant N to A weak 
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no. confs 


K or ami 


am % 


aM N 


aM A 


Quenched 32 3 x 64 a" 1 = 2. 14(6) GeV 


200 
200 
200 


0.1554 
0.1558 
0.1562 
k c =0.1571 


0.263(2) 
0.229(2) 
0.192(2) 
0. 


0.592(5) 
0.556(6) 
0.518(6) 
0.439(4) 


0.687(7) 
0.666(8) 
0.646(9) 
0.598(6) 


Unquenched Wilson 24 3 x 40 [@] a~ l = 2.56(10) GeV 


185 
157 


0.1575 
0.1580 


0.270(3) 
0.199(3) 


0.580(7) 
0.500(10) 


0.645(5) 
0.581(14) 


Unquenched Wilson 24 3 x 32 [§] a" 1 = 2.56(10) GeV 


200 


0.15825 
Kc = 0.1585 


0.150(3) 
0. 


0.423(7) 
0.366(13) 


0.533(8) 
0.486(14) 


MILC 20 3 x 64 a -1 = 1.58 GeV 


150 
150 


0.03 
0.02 


0.373(3) 
0.306(3) 


0.886(7) 
0.800(10) 


1.057(14) 
0.992(16) 


MILC 28 3 x 64 cT x = 1.58 GeV 


118 


0.01 


0.230(3) 


0.751(7) 


0.988(26) 



Table 1: The number of configurations, the hopping parameter, fc, for the case of Wilson fermions or the 
mass of the light quarks, m/, for the case of staggered quarks, the pion, nucleon and A mass in lattice units. 



matrix element can be expressed in terms of four transition form factors as 



1/2 



u x {p',s) 



{gliigpv-gXpguv) q P +Ci{q 2 )gxp + 6 2 a X a fi 



M 2 N 



u(p,s) (3.1) 



where = p' ^ —Pfi is the momentum transfer and A 3 (x) = \\f{x)y^y^\if{x) is the isovector part 
of the axial current (t 3 being the third Pauli matrix). In order to evaluate this matrix element on the 
lattice we compute the three point function (G^ N (t2,h\p' ; P;T)) . We eliminate the exponential 
decay in time and the overlaps of the interpolating fields with the physical states by forming an 
appropriate ratio, R a (t2,ti;p' ,p;T; jj.), of three-point and two-point functions given by 



Re 



<G£(fe,p';r 4 )) 



(G w (f 2 - ?1 ,p;r 4 )) <Gg(ti,p , ;r 4 )> (G£(* 2 y ; r 4 )> 



(Gf^-tup 1 ;^)) (G^(fi,p;r 4 )> (G N (t 2 ,p;r 4 )) 



1/2 



n ff (p / ,p;P,M) 



(3.2) 



where T 4 




and r , 



Oj0 




With t\ we denote the time when a photon interacts 



with a quark and with ?2, the time when the A is annihilated. The ratio given in Eq. (3.2) is 
constructed so that the two-point functions that enter have the shortest possible time separation 
between source and sink. This provides an optimal signal to noise ratio. For large time separations 
t\ and t2, the ratio R a (t2,h;p' ,p;T; /i) becomes time independent and yields the transition matrix 
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element of Eq. ( fO| ) up to the renormalization constant Za- The latter has been computed non- 



perturbatively using the RI-MOM method for quenched Q13Q and two flavor of dynamical Wilson 
fermions [14]. The values obtained in both cases are all consistent with Za = 0.8. For domain wall 
fermions we use the values given in Ref. [j|]. We use kinematics where the A is produced at rest 
and by Q 2 = —q 2 we denote the Euclidean momentum transfer squared. 

There are various choices for the Rarita-Schwinger spinor index a and projection matrices T 
that yield the four axial form factors. Each of these choices requires a separate sequential inversion. 



As in the case of the evaluation of the electromagnetic N to A transition form factors [11] we use 
optimized A sources in order to maximize the number of lattice momentum vectors contributing to 
a given Q 2 value. The optimized A sources turn out to be the same as those used in our study of the 
electromagnetic form factors [fPT[]. Namely we use the combinations Si (q; jj.) = £a=i n^^T/i; ju), 



52(q;/i)=I^ fc=1 n ff (q;r fc ;/x)and5 3 (q;/x) = n3(q;r 3 ;ju)-I(ni(q;ri;ju) + n 2 (q;r 2 ;/x)). The 
four axial form factors can be extracted from the following expressions 



5i(q;j) = iB 



2(E N + M N ) 



| (E N + M N ) (E N - 2M A + Mff) + ( £ pt\ p j | 



MN ^-M A )Ci + M N Ci-^ NP j(i/ 



7 = 1,2,3 



5!(q;4) =fl£p* 

k=\ 



r A,MA A 



En-Ma^a 



S 2 (q;j) = iA 
5 3 (q;j) = iA 



I P k {S fl (p 2 - P 3 ) + 8 j2 (p 3 -p i ) + 8 p {p l - p 2 )) C 



-[S n (p 2 p 3 )-8 j2 (p l p 3 ))cj 



(3.3) 



where A = ^/2/3 yjE N /(E N + M N )/( 3E N M N ) and B =A/(E N +M N ). The axial form factors can 



be extracted by performing an overconstrained analysis as described in Refs. [11, 12]. 



4. Results 

In Fig. [T] we show our lattice results for the four axial form factors for quenched and un- 
quenched Wilson fermions and in the hybrid approach. We observe that C3 is consistent with 
zero and that unquenching effects are small for the dominant form factors, C5 and Cg . The form 
factor C4 shows an interesting behavior: The unquenched results for both dynamical Wilson and 
domain wall fermions show an increase at low momentum transfers. Such large deviations between 
quenched and full QCD results for these relatively heavy quark masses are unusual making this an 
interesting quantity to study effects of unquenching. 

In the chiral limit, axial current conservation leads to the relation C^(Q 2 ) = M^C^ (Q 2 ) / Q 2 . In 
Fig. U we show the ratio (Q 2 /M^) C^(Q 2 )/C^(Q 2 ) for quenched and unquenched Wilson fermions, 
and in the hybrid scheme. In each case we show results for the available quark masses and in the 
chiral limit. The expected value in the chiral limit for this ratio is one. For finite quark mass the 
axial current is not conserved and for Wilson fermions chiral symmetry is broken so that deviations 
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Figure 1: The axial form factors C3 , Cf , Cg and Cg as a function of Q 2 . In all plots we show quenched re- 
sults, denoted by Nf = 0, at K = 0.1554 (crosses), at K = 0.1558 (open circles) and at K = 0.1562 (asterisks). 
The graphs on the left hand side also show unquenched Wilson results, denoted by Nf = 2, at K = 0.1575 
(triangles), K = 1580 (filled circles) and K = 0.15825 (open squares). The graphs on the right hand side also 
show results in the hybrid approach at am/ = 0.03 (crosses), ami — 0.02 (filled diamonds) and ami = 0.01 
(filled triangles). 

from one are expected. We observe that this ratio differs from unity at low Q 2 but approaches unity 
at higher values of Q 2 . For the hybrid scheme the ratio is consistent with unity even at the lowest 
available Q 2 , as expected for chiral fermions. That such chiral restoration is seen on the lattice even 
when using Wilson fermions demonstrates that lattice methodology correctly encodes continuum 
physics. 

At finite pion mass partial conservation of axial current (d^A" (x) = f K m 2 n a (x)) leads to the 
non-diagonal Goldberger-Treiman relation C^(Q 2 ) = fngxNAiQ 2 ) /2Mn where gxNA(Q 2 ) is deter- 
mined from the matrix element of the pseudoscalar density 

2m q <A + (p',s')\xj/Y5^-V\N(p,s) >= g KNA (Q 2 )\^- q ° 2 »a{p' ,s')u{p,s), (4.1) 

2 V 3 2M N Q z + m A n 

where m q is the renormalized quark mass. The pion decay constant, f K , is determined from the 
two-point function < 0\A a ^\7i h {p) >= ip^8 ah f n , defined so that the continuum value is taken 
f K = 93.2 MeV. In order to relate the lattice pion matrix element to its physical value we need 
the pseudoscalar renormalization constant, Z p . We take for quenched [ |l3| ] and dynamical Wilson 
fermions jl4| ] Z p (/x 2 a 2 ~ 1) = 0.5(1) computed using the RI-MOM method. This value may de- 
pend on the renormalization scale whereas it is not known for domain wall fermions. In Fig. || we 
show the result for g K NA for Wilson fermions and the linear extrapolation in m 2 n of these results to 
the chiral limit. We note that in the figure we only show statistical errors which do not include a 
10% uncertainty in Z p . Furthermore, we would like to stress that the determination of the quark 
mass m q using the axial Ward identity has corrections of order a. These corrections become more 
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Figure 2: The ratio (Q 2 /Ml)C^/C^ versus Q 2 . 
Top: In the quenched theory, Middle: For dynam- 
ical Wilson fermions, Bottom: In the hybrid ap- 
proach. 
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Figure 4: The ratio fng K NA/ {2M N C^) ver- 
sus Q 2 for quenched and dynamical Wilson 
fermions. 
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Figure 3: g K NA as a function of Q 2 for quenched 
and dynamical Wilson fermions. Top: For our 
three- K values. Bottom: In the chiral limit. 
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Figure 5: The ratio C^jC\ as a function of Q 2 for 
quenched QCD, for dynamical Wilson fermions and 
for the hybrid scheme. 



significant with decreasing quark mass. These can lead to large uncertainties in the determination 
of g TIN A- 

In Fig. U we show the ratio fngnNA/ (2M N C^) for Wilson fermions. As can be seen this ratio 
is almost Q 2 independent and as the quark mass decreases it becomes consistent with unity in 
agreement with the non-diagonal Goldberger-Treiman relation. 

Under the assumptions that C3 ~ and that C4 is suppressed as compared to Cf, both of which 
are justified by the lattice results, the parity violation asymmetry can be shown to be proportional 



to the ratio C^/C 3 V [0. The form factor can be obtained from the electromagnetic N to A 
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transition. Using our lattice results for the dipole and electric quadrupole Sachs factors, 54/1 and 
&E2 0, C 3 y is extracted from the relation C\ = f (&mi -^ei). 

We show in Fig. |5| the ratio /C\ , for pions of mass about 500 MeV. As can be seen un- 
quenching effects are small and in order to assess the quark mass dependence we extrapolate our 
quenched results, which carry the smallest errors, to the chiral limit. We find only a small increase 
in this ratio as we tune the quark mass to zero, indicating a weak quark mass dependence. There- 
fore our lattice evaluation provides a prediction for the physical value of this ratio, which is the 
analog of the gA/gv- Our lattice results show that this ratio, and therefore to a first approximation 
the parity violating asymmetry, is non-zero at Q 2 = and increases for Q 2 ~ 1.5 GeV 2 . 

5. Conclusions 

In summary we have provided a lattice calculation of the axial N to A transition form factors in 
the quenched approximation, using two degenerate dynamical Wilson fermions and within a hybrid 
approach where we use MILC configurations and domain wall fermions. 

The main conclusions are: 1. is consistent with zero whereas C4 is small but shows the 
largest sensitivity to unquenching effects. 2. The two dominant form factors are C5 and Cg. 
These are related in the chiral limit by axial current conservation. The ratio (Q 2 /M^C^/C^, 
which must be unity if chiral symmetry is unbroken, is shown to approach unity as the quark mass 
decreases. 3. For any quark mass the strong coupling g K NA and the axial form factor Cj show 
a similar Q 2 dependence with the non-diagonal Goldberger-Treiman relation being reproduced as 
the quark mass decreases. 4. The ratio of C5 /C^ which determines to a good approximation the 
parity violating asymmetry is predicted to be non-zero at Q = and has a two-fold increase when 
Q 2 ~ 1.5 GeV. 
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